Relativistic models of the neutron-star matter equation of state 
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Motivated by a recent astro-physical measurement of the pressure of cold matter above nuclear- 
matter saturation density [l] , we compute the equation of state of neutron star matter using accu- 
rately calibrated relativistic models. The uniform stellar core is assumed to consist of nucleons and 
leptons in beta equilibrium; no exotic degrees of freedom are included. We found the predictions of 
these models to be in fairly good agreement with the measured equation of state. Yet the Mass-vs- 
Radius relations predicted by these same models display radii that are consistently larger than the 
observations. 
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The quest for the Holy Grail of Nuclear Physics — the 
equation of state (EOS) of hadronic matter — remains an 
area of intense activity that cuts across a variety of dis- 
ciplines. Indeed, the limits of nuclear existence, the dy- 
namics of heavy- ion collisions, the structure of neutron 
stars, and the collapse of massive stellar cores all de- 
pend sensitively on the equation of state. With the ad- 
vent and commissioning of sophisticated new radioactive 
beam facilities, powerful heavy-ion colliders, telescopes 
operating at a variety of wavelengths, and more sensitive 
gravitational wave detectors, one will be able to probe 
the nuclear dynamics over a wide range of nucleon asym- 
metries, temperatures, and densities. However, in the 
present contribution we focus on the dynamics of cold 
matter under extreme conditions of density (both small 
and large) and for this case neutron stars remain the tool 
of choice - Being both very compact and extremely 
dense, neutron stars are unique laboratories for probing 
the equation of state of neutron-rich matter under con- 
ditions unattainable by terrestrial experiments. 

Intimately connected to the equation of state of cold, 
neutron-rich matter is the mass-vs-radius (M-R) relation- 
ship of neutron stars. Indeed, an EOS is the sole ingre- 
dient that must be supplied to solve the equations of 
stellar structure (i.e., the Tolman-Oppcnhcimcr-Volkoff 
equations). Conversely, knowledge of the M-R relation is 
sufficient to uniquely determine the equation of state of 
neutron star matter [5]. As argued by Lindblom almost 
20 years ago, the availability of such information — even 
from a single neutron star — will provide interesting in- 
formation about the equation of state [5] . Viewed in this 
light, the recent report of combined mass-radius mea- 
surement for three neutron stars and the subsequent de- 
termination of the equation of state is significant [T] . In 
particular, the conclusion that the EOS so determined 
is softer than those containing only nucleonic degrees of 
freedom is both interesting and provocative. 

In this contribution we compute the equation of state 
of neutron star matter and the resulting M-R relation 



using accurately-calibrated relativistic mean-field mod- 
els. These models have been calibrated to the properties 
of infinite nuclear matter at saturation density [5J, to 
the ground-state properties of finite nuclei [7j |8] , or to 
both [9]. Unlike the former two, the latter parametriza- 
tion predicts a significantly soft symmetry energy, a fea- 
ture that appears consistent with the behavior of dilute 
neutron matter (see Ref. [TU] and references therein). A 
detailed explanation of the role of the model parameters 
on the equation of state is given below. We note, how- 
ever, that none of the models considered in this work in- 
clude exotic degrees of freedom, such as hyperons, meson 
(condensates), or quarks. In this regard, our results are 
mixed when compared with the conclusions of Ref. pQ. 
On the one hand, the stellar radii predicted by the rela- 
tivistic models are larger than observed, seemingly con- 
firming that such equations of state are too stiff. On the 
other hand, the agreement between the predicted and 
observed EOS suggests the opposite. 

The structure of neutron stars is sensitive to the equa- 
tion of state of cold, fully catalyzed, neutron-rich mat- 
ter over an enormous range of densities [2J-0]. For the 
low-density outer crust we employ the equation of state 
of Baym, Pcthick, and Sutherland [TT]. At densities of 
about a third to a half of nuclear-matter saturation den- 
sity, uniformity in the system is restored and for this 
(liquid-core) region we use an EOS derived from a rep- 
resentative set of accurately calibrated relativistic mean- 
field models [BHS]- It nas been speculated that the re- 
gion between the outer crust and the liquid core con- 
sists of complex and exotic structures, collectively known 
as nuclear pasta [I2TH4] . Whereas significant progress 
has been made in simulating this exotic region [15H17) . 
a detailed equation of state is still missing. Hence, we 
resort to a fairly accurate polytropic EOS to interpo- 
late between the solid crust and the uniform liquid inte- 
rior [TU [19]. To compute the transition density from 
the liquid core to the solid crust we employ a rela- 
tivistic random-phase-approximation (RPA) analysis to 
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TABLE I: Parameter sets for the four models used in the text to generate the equation of state. The parameter n and the 
meson masses m s , m v , and m p are all given in MeV. The nucleon mass has been fixed at M = 939 MeV in all the models. 



search for the critical density at which the uniform sys- 
tem becomes unstable to small amplitude density oscil- 
lations |19j . 

Accounting for most of the stellar radius and practi- 
cally all of its mass, the liquid core is structurally the 
most important component of the star. Matter in the 
liquid core is assumed to be composed of neutrons, pro- 
tons, electrons, and muons in chemical equilibrium. We 
reiterate that no exotic degrees of freedom are included 
in the model. Both electrons and muons are treated as 
non-interacting relativistic Fermi gases. For the hadronic 
component, the equation of state is generated using 
accurately-calibrated relativistic models. Details on the 
calibration procedure may be found in Refs. [50Tl23| . 
The model includes a nucleon field (ip) interacting via 
standard Yukawa couplings to two isoscalar mesons (a 
scalar (j) and a vector V 1 ) and one vector-isovector me- 
son (6 M ) £321 [21] . Such an interacting Lagrangian density 
may be written as follows [5J [50] [5T] : 
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In addition to the Yukawa couplings (g s , g v , and g p ), the 
model is supplemented by non-linear meson interactions 
given by 
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The inclusion of scalar cubic (n) and quartic (A) self- 
interactions dates back to the late seventies [5J] and is in- 
strumental for softening the incompressibility coefficient 
of symmetric nuclear matter, as required to explain the 
excitation of the the nuclear breathing mode [2 5) . 

Of particular interest and of critical importance to the 
present study are the vector self-interaction (£) and the 
isoscalar-isovector mixing term A v [pj [55] . That both of 
these parameters are zero in the enormously successful 
NL3 model (see Table [I]) suggests that existing laboratory 
data are fairly insensitive to the physics encoded in these 
two parameters. Indeed, Muller and Serot found possible 



to build models with different values of £ that reproduce 
the same observed properties at normal nuclear densities, 
yet produced maximum neutron star masses that differ 
by almost one solar mass [BJ. This result indicates that 
observational data on neutron stars — rather than labora- 
tory experiments — may provide the only meaningful con- 
straint on the high-density component of the equation of 
state. Further, it indicates that the empirical parameter 
C provides an efficient tool to control the high-density 
component of the equation of state. 

The isoscalar-isovector coupling constant A v was intro- 
duced in Ref. [55] to modify the poorly known density de- 
pendence of the symmetry energy. The symmetry energy 
represents the energy cost involved in changing protons 
into neutrons (and vice- versa). To a good approxima- 
tion, it is given by the difference in energy between pure 
neutron matter and symmetric nuclear matter. With 
only one isovector parameter (g p ) to adjust, relativis- 
tic mean-field models have traditionally predicted a stiff 
symmetry energy. The addition of A v provides a simple — 
yet efficient and reliable — method of softening the sym- 
metry energy without compromising the success of the 
model in reproducing well determined ground-state ob- 
servables [9]. Indeed, whereas models with different val- 
ues of A v reproduce the same exact properties of symmet- 
ric nuclear matter, they yield vastly different predictions 
for both the neutron radii of heavy nuclei and for the 
radius of neutron stars [55] [5BJ . Given that the neutron 
star radius is believed to be primarily controlled by the 
symmetry pressure at intermediate densities [4], the up- 
coming Parity Radius Experiment (PREx) at the Jeffer- 
son Laboratory (with an imminent start date of March, 
2010) will provide a unique laboratory constraint on a 
fundamental neutron star property |27l I28j . 

In summary, the two empirical parameters £ and A v 
provide a highly economical and efficient control of the 
softness of the high-density component of equation of 
state and of the symmetry pressure at intermediate densi- 
ties, respectively — with the former primarily controlling 
the maximum neutron star mass and the latter the stel- 
lar radius. Parameter sets for all the models employed in 
this work are listed in Table H 
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FIG. 1: (Color online) Mass-vs-Radius relation predicted by 
the four relativistic mean-field models discussed in the text. 
The observational data represent la confidence contours for 
the three neutron stars reported in Ref. [I] 



In Fig. [T] we compare observational results for three 
neutron star masses and radii against the model predic- 
tions. These neutron stars are in the binaries 4U 1608- 
52 [29], EXO 1745-248 30], and 4U 1820-30 [31]. The 
very stiff behavior of the NL3 equation of state is imme- 
diately evident. With both empirical parameters C an d 
A v set equal to zero, it is not surprising that the NL3 
model predicts neutron star masses as large as 2.8 Mq 
with very large radii. As compared to the observational 
data, the NL3 model suggests a radius for a 1.7 solar- 
mass neutron star that is about 6 km too large. More- 
over, the NL3 equation of state is so stiff that gravity 
in a 2.8 Mq neutron star can compress matter to only 
about four times normal nuclear density (see Table |ll| . 
All these, even when the model provides an excellent de- 
scription of many laboratory observables. 

As first suggested by Muller and Serot [BJ, adding a 
vector self-interaction (with £ = 0.06) dramatically re- 
duces the repulsion at high densities and ultimately the 
limiting neutron star mass. As compared to the NL3 pa- 
rameter set, the maximum neutron star mass predicted 
by Muller and Serot (MS) is reduced by almost one solar 
mass (see Fig. [IJ and Table [IIJ) . Consistent with this soft- 
ening is a significant increase in the compactness of the 
star. For example, for a neutron star mass of 1.8 Mq, 
NL3 predicts a stellar radius that is more than 3 kilo- 
meters larger than MS. Note, however, that the density 
dependence of the symmetry energy predicted by NL3 
and MS is practically identical (see inset in Fig. [2]) . In 
particular, this is reflected in the identical prediction of 
0.28 fm for the neutron-skin thickness of 208 Pb. This 
suggests that tuning the density dependence of the sym- 
metry energy — via the addition of the isoscalar-isovector 
mixing term A v — may yield a further reduction in neu- 
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TABLE II: Predictions for the central baryon density, central 
pressure, mass, and radius of the limiting neutron star for the 
four models employed in the text. The last column lists pre- 
dictions for the radius of a 1.4 solar-mass neutron star. The 
baryon density is given in fm -3 , the pressure in MeV fm -3 , 
the mass in solar masses, and the radii in kilometers. 



tron star radii [26], as suggested by observation. 

Incorporating information on nuclear collective modes 
in the calibration procedure of the FSUGold model favors 
a non-zero value for A v [5]. Further, it now seems that 
the resulting softening of the symmetry energy is consis- 
tent with the EOS of dilute neutron matter predicted by 
various microscopic approaches (see Refs. [T0"H3"2"H3~4"] and 
references therein) . That the addition of A v produces the 
intended effect can be appreciated in Fig. [I] and Table |H] 
That is, although one has adopted the same value of C for 
both MS and FSUGold, their predictions for the radius of 
a "canonical" 1.4 solar-mass neutron star differ by more 
than one kilometer. Related to this fact is the signifi- 
cantly smaller neutron-skin thickness of 208 Pb predicted 
by FSUGold (0.21 fm vs 0.28 fm). However, it appears 
that the combined softening of the EOS at high densities 
(through £) and of the symmetry pressure (through A v ) 
is insufficient to explain the observational data; the min- 
imum stellar radius predicted by the FSUGold model is 
about 11 km, significantly larger than suggested by ob- 
servation. 

In an effort to describe the observational data, we have 
constructed an "Extra Soft" (XS) relativistic mean-field 
model constrained by the properties of symmetric nu- 
clear matter at saturation density (i.e., equilibrium den- 
sity, binding energy per nucleon, and incompressibility 
coefficient). In regards to these properties, the model 
is indistinguishable from FSUGold. The only additional 
constraint imposed on the model is that its limiting mass 
be no smaller than 1.6 solar masses. We feel that lower- 
ing this limiting value any further may start conflicting 
with the observational data. Although no exhaustive pa- 
rameter search was conducted, we trust that the result- 
ing extra-soft equation of state (as given in Table [1} is 
representative of the softness that may be achieved with 
present-day relativistic mean-field models. With such a 
soft model, neutron star radii get significantly reduced 
indeed (see Fig. [IJ). For example, the neutron radius of 
a 1.4 Mq neutron star is reduced by almost one kilo- 
meter relative to the FSUGold prediction (see Table [TTJ) 
and by more than 1.5 km at its limiting mass of 1.6 Mq. 
Still, the minimum neutron star radius of R = 10.41 km 
predicted by the model remains outside the reported la 
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FIG. 2: (Color online) Equation of state (Pressure vs baryon 
density) of neutron star matter predicted by the four rela- 
tivistic mean-field models discussed in the text. The three 
data point are from the observational extraction as described 
in Ref. [JJ. The symbols (stars) indicate the central density 
and pressure for the maximum-mass neutron star. The inset 
shows the symmetry pressure, given as the pressure of pure 
neutron matter minus that of symmetric nuclear matter. 

confidence contours pQ. 

Do we then conclude that the results presented in 
Fig. [I] are indicative of relativistic equations of state that 
are too stiff? Do the observational results unambigu- 
ously called for a softer equation of state, as would be 
produced by exotic states of matter, such as meson con- 
densates and/or quark matter? To answer this question 
we compare in Fig. [2] the various equations of state used 
to generate Fig. [I] against the values extracted from the 
observational data pp. The inset in the figure displays 
the symmetry pressure for the models under considera- 
tion. The observed softening of the symmetry pressure 
between models is entirely due A v . Note, however, that 
unlike the neutron-skin thickness of neutron-rich nuclei, 
the radius of a neutron star is not uniquely constrained 
by the symmetry pressure at low to intermediate densi- 
ties |26j . Thus models with similar symmetry pressures 
may — and do — predict significantly different stellar radii. 
Contrary to the expectations generated by Fig. [I] most 
of the equations of state are not too stiff. Indeed, with 
the exception of NL3, the remaining equations of state 
appear, if anything, slightly too soft at the highest den- 
sity. Based on these results — and these results alone — 
nucleonic equations of state do not seem to be in conflict 
with the observational data. 

In summary, the Mass-vs-Radius relation of neutron 
stars was computed using equations of state derived from 
relativistic mean-field models. Although the models are 
calibrated in the vicinity of nuclear-matter saturation 
density, it is possible to tune their high-density behav- 
ior in a highly efficient and economical manner. In this 
contribution we have used two parameters to control the 



maximum neutron star mass and the stellar radius. As 
we compared our predictions to the observational data 
a conflict emerged. Whereas one could generate equa- 
tions of state that are in agreement with observation, the 
predicted stellar radii are too large. This result is par- 
ticularly intriguing given that "inversion" methods exist 
for extracting the equation of state of stellar matter di- 
rectly from masses and radii of neutron stars (Sj. Thus, 
one would expect that if the M-R predictions do not 
match observation, neither would the equations of state. 
Clearly, to reconcile these facts much work remains to be 
done in both the observational and theoretical fronts. For 
now we must conclude — although the existence of exotic 
stars is very appealing — that the downfall of the purely 
nucleonic equations of state may be premature. 
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